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Abstract 

We present numerical evidence for the hypothesis that, in the planar limit, four dimensional Eu- 
clidean Yang-Mills theory on a finite symmetrical four-torus breaks chiral symmetry spontaneously 
when the length of the sides I is larger than a critical value l c with a bilinear condensate whose 
value is independent of I. Therefore spontaneous symmetry breaking occurs at finite volume and 
infinite N P reduction holds for the chiral condensate. 
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I. INTRODUCTION. 



The sum of pure gauge connected vacuum Feynman diagrams goes as iV c 2 at fixed 't Hooft 
coupling, A = g 2 N c [lj. Fermions make a contribution that goes only as N c . Similar behavior 
occurs in the strong coupling expansion on the lattice. In short, the mere fact that there are 
order N 2 gluons and only order N c fermions indicates that at leading order in N c fermions 
propagate in a medium determined by pure gauge theory, with no "back reaction" from the 
fermions. 

Until recently, QCD oriented computer simulations were done in the so called quenched 
(valence) approximation where the fermionic determinant is set to unity. Thus, the fermionic 
"back reaction" is eliminated by hand. Results, perhaps surprisingly, came out quite close to 
experiment. One way to understand this is to argue as in the previous paragraph, and then 
assert that this specific property of the planar limit holds approximately also at N c = 3. 
The valence approximation is not self-consistent for any finite N c and does not make up 
an acceptable field theory. The main problems occur for massless quarks: the fermion de- 
terminant can be zero exactly for massless quarks and there are symmetries obeyed by the 
massless fermion propagator but violated by the determinant (anomalies). These properties 
are important and get lost when the fermion determinant is replaced by unity. If in Na- 
ture all quarks were very massive relative to Aqcd ; the quantitative success of the valence 
approximation would be easy to understand without appealing to the large N c expansion. 
But, the up and down quarks are much lighter than Aqcd and dominate the low energy 
physics of strong interactions. 

Configurations for which the fermion determinant is exactly zero can be eliminated, thus 
redefining the quenched approximation to mean that the fermion determinant is unity for 
gauge fields of zero topology and zero for all other gauge fields. This augmented quenched 
approximation can be easily implemented on the lattice if one uses exactly chiral lattice 
fermions. Furthermore, simulations in the large N c limit with massless quarks in the fun- 
damental representation should be restricted to zero topology: The fermion determinant is 
zero for non-zero topology and then it obviously overcomes the pure gauge action exponent, 
but it is much closer to unity than the exponent of the pure gauge action for zero topology 
and therefore can in this case be replaced by unity at leading order in For practical 
simulations the lengths of times are such that lattice gauge configurations in one topological 
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sector do not evolve into another sector when N c is large enough (depending on the lattice 
volume, N c larger than order 10 eliminates topology changes). Thus, even for massive quarks 
one ends up working always at a unique total topology. The other problem one has with 
the quenched approximation at finite N c , that of anomalies, becomes explicitly an effect of 
subleading order in the expansion. At finite N c the inconsistency of the quenched approx- 
imation stemming from the absence of the anomaly is reflected by new chiral divergences. 
These divergences are absent when fermion loops are included. But these divergent terms 
have coefficients that explicitly vanish in the large N c limit, validating the self-consistency of 
the leading term in the large N c expansion. In short, at N c = oo, fermions are automatically 
"quenched", and the field theoretical framework is intact. 

In the context of QCD quenching was introduced for the simple reason that computers 
were not powerful enough to include the fermion determinant, while excluding this determi- 
nant produced a model that one could simulate in practice. One would guess therefore that 
relatively modest computer resources (by the standards of today) would suffice to numer- 
ically solve for, say, the meson spectrum of large N c gauge theory. However, going to the 
large N c limit increases computational cost, so additional insights into simplifications that 
occur at infinite N c are needed to find ways to offset this. 

One simplification is that one does not have to worry about attaining the "thermodynamic 
limit"- in the usual sense of the word - when N c is infinite : When N c = oo any Wilson 
loop expectation value on a finite Z 4 torus does not depend on I, so long as I > l c , where l c is 
a fixed physical length, of order | fermi in QCD terms. Another simplification is afforded by 
a specific trick, first introduced by Gross and Kitazawa jj] and later adapted to the lattice 
and tested in two Euclidean dimensions p , which eliminates finite I effects also from meson 
propagators. The essence of this trick is to "force-feed" a differential of momenta into the 
two quark lines that make up the meson propagator. 

A putative puzzle now presents itself threatening to invalidate the exact independence of 
meson propagators on I. Obviously, these propagators are strongly affected by spontaneous 
chiral symmetry breaking for massless quarks at any N c , including N c = oo. But, symmetries 
tend not to break spontaneously in a finite volume, as can easily be checked explicitly by 
solving A($ 2 ) 2 at infinite N, where N is the number of components of the real scalar field $ 
and the calculation is done on a finite four torus [6] . The resolution of the puzzle is simple: 
unlike in the case just mentioned, spontaneous chiral symmetry breakdown is induced by 
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the N c — > oo limit even on a finite torus, of side /, at least so long as I > l c . Moreover, 
the value of the traditional order parameter, the expectation value of the fermion bilinear 
condensate, is / independent in any well defined renormalization scheme that in itself is 
volume independent. 

This resolution of the puzzle is the central hypothesis of the present paper. We shall 
present numerical evidence verifying it. Also, we shall argue that spontaneous symmetry 
breaking is very natural at infinite N c , strengthening insights developed during the last 
decade on the connection of chiral symmetry breaking with random matrix theory. In short, 
spontaneous chiral symmetry breaking occurs at infinite N c as a result of sufficient disorder, 
which seems to always be present so long as we have confinement, at I > l c , but certainly 
does not require confinement. 

II. SPONTANEOUS CHIRAL SYMMETRY BREAKING ON THE LATTICE. 

Learning how to exactly preserve chiral symmetry on the lattice has allowed to separate 
the understanding of the phenomenon of spontaneous chiral symmetry breaking (S^SB) 
into two questions. The first asks whether S^SB occurs at fixed lattice spacing and by 
what mechanism. The second asks whether the effect survives the process of taking the 
continuum limit, where the lattice spacing goes to zero. Because ultraviolet physics seems 
largely irrelevant for the infrared phenomenon of S^SB it is widely believed that the first 
question addresses the more fundamental dynamical issues. 

We start our study focusing on the first question. At this stage our simulations are carried 
out at one - relatively coarse - lattice spacing a, which, in QCD terms is about | fermi 
and at zero total topology. Addressing the first question is the main objective of this paper. 
Subsequently, we proceed to the second question and check scaling, which is our secondary 
objective here. The second question is not fully addressed because we do not determine the 
needed mass normalization nonperturbatively. In addition, the amount of data collected at 
smaller lattice spacings is relatively small. 

Our quarks are massless and we should only consider zero topology. However, the tech- 
nique we use to establish S^SB extends to a prediction at non-zero topology, providing a 
crosscheck. We crosscheck our results against a few simulations at nonzero topology. Also, 
since we really define the large N c limit by taking N c to infinity before taking the quark 
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mass to zero, the restriction to zero topology might be questioned. We check whether our 
zero quark mass results for the chiral condensate are compatible with those at finite quark 
mass. 

In order to get some feeling whether the numerical behavior that we see is reasonable we 
compare our methods to the situation in two dimensions, the 't Hooft model. There the 
exact value of the condensate in the continuum limit is known, so we can independently 
judge how far astray the simulations might lead us. 

The bilinear fermionic action is described by the massless overlap Dirac operator []| D : 



D„ 



i+v 
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V- 1 = = 75^75 = sign(^(M)) 75 (1) 

H W (M) is the Wilson Dirac operator at mass M, which we shall choose as M = — 1.5. M 
should not be confused with the bare quark mass m, to be introduced below. 



H W (M) = 75 



M + 4 - Y ( ^—^T, + ^^Tl 



(2) 



The Tjj matrices are the lattice generators of parallel transport and depend parametrically 
and analytically on the lattice links U^(x) which are SU(N C ) matrices at site x associated 
with the link connecting site x to site x + p,, where fi is a unit vector in the positive \i- 
direction. 

At finite N c and with Nf flavors we would have to take into account a factor of [det D a ] f 
influencing the distribution of the gauge fields. As explained in the introduction, at zero 
topology and infinite N c this factor is replaced by unity. The internal fermion-line propaga- 
tor, —^7 is therefore not needed. For fermion lines continuing external fermion sources we 

nn 

are allowed to use a slightly different quark propagator 8, 9J defined by: 

A ~ 1 + V [6) 

A = —A^ and anticommutes with 75. The spectrum of A is unbounded, but is determined 
by the spectrum of V which is restricted to the unit circle. One should think of A as 
dimensionless, and of \M\ as providing the needed dimension. Up to a dimensionful unit, A 
should be thought of as a lattice realization of the continuum massless Dirac operator, D: 

2\M\A ^D = 7^ + (4) 
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A positive, dimensionless, quark mass is added by defining 



A(jm)=A + h (5) 

The relation between the bare, properly normalized, quark mass m and the parameter /i is 
m = 2|M|/i. This can be used to define a massive overlap Dirac operator, D Q (fi), appropriate 
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for the internal fermion lines making up the fermion determinant 

The eigenvalues of V can be written as — e~ 2ld . Then the eigenvalues of A are ztan#. 
The eigenvalues of the massless overlap operator D are given by ze sin0. For \9\ « 1 
the eigenvalues of A and D are numerically very close. 

At finite N c the occurrence of S^SB manifests itself mathematically as a lack of commu- 
tativity of the limits V — > oo (here V is the lattice volume) and fi —>■ 0. 

lim^o limbec y{TrA- 1 (fj l )) NcjV = E/ 

lim^ (TrA-\fi)) Nc y = (6) 

(...) means gauge averaging, which includes the factor [det D (n)] N f] and this is an additional 
source of dependence on /i. The subscript "iV c , V" means that the average is performed on 
the torus of volume V at fixed and finite N c . 

As N c — > oo, S will diverge, but E = jj- will have a finite limit. The hypothesis of this 
paper is that: 

lirn^o hmAr c -*oo ^-(TrA _1 (^))jv c ,v = S ^ 

Um^ (TrA- l {fi)) Nc y = (7) 

The second line of the above two sets of equations is a direct consequence of 75A + A75 = 
0. The notation also implies that S is independent of the volume V, meaning that 
lim^o h m A r c -+oo jf(TrA~ l (n))]s[ c y has to come out a linear function of V with zero in- 
tercept. A somewhat superficial (but nevertheless quite convincing) argument supporting 
this claim starts by imagining an expansion of ■^-(TrA~ 1 (fi))N c ,v i n traces of Wilson loops, 
and then invokes large N c reduction for the latter. However, the argument ignores the ques- 
tion whether the expansion in loops converges - a substantial issue in particular since we 
intend to take the bare quark mass to zero at the end. 
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'he quantity X, as defined above, diverges in the quenched approximation for finite 
1(3, [llj]. However, the quantity S, as defined above, stays finite in the N c — > 00 limit. 



We learn that we need to refine our definition of the large N c limit for the case of massless 
quarks. One needs to take the limits in the following order: First one takes N c — *■ oo in 
the 't Hooft way, at fixed number of flavors, Nf > 1, and a fixed mass, /i ^ 0, and only 
subsequently one lets \i — *■ 0. The first step removes the dependence on Nf. There exist 
other orderings of limits, in which the dependence on Nf can be preserved. With our order 
of limits, the independence on Nf means that one could obtain the same limit in yet another 
way: First one takes N c —>■ oo in the 't Hooft way, for the quenched theory at a fixed mass 
H > 0, and only subsequently one lets // — > 0. This is the basic definition of the large N c 
limit we adopt in this paper. The spectra of the massless and massive Dirac operator are 
trivially related in the quenched theory. Therefore, we shall see later that if we are careful 
about what we are doing we can get away by taking the large N c limit of the quenched 
theory at fi = directly; what exactly we do and why will become clear only after we speak 
in greater detail about infrared divergences coming from quenching. 



III. RANDOM MATRIX THEORY (RMT) AND S^SB 

Clearly, the noncommutativity of limits described in the previous section implies nontriv- 
ial crossover effects in the dependence of the quantity (TrA~ 1 (fi)) Nc y on fi, V and N c . This 
rapidly becomes a numerical burden, since isolating the right regime might require extremely 



large values of V_ 
many workers 



fat fixed N c ) or N c (at fixed V). Fortunately, due to the successful efforts of 
12J, the physics of the crossover regime is quite well understood. As is usually 
the case, the crossover is largely universal, and, in this case, there exists a representative of 
its universality class given by a Gaussian random matrix model, introduced by Shuryak and 
Verbaarschot jlHI ]. 

For Nf > 1 and finite N c the crossover regime controlling the pair of limits // — > and 
V — > oo can be described by traditional methods, employing chiral Lagrangians jl^j . At 
zero topology, the dependence on fi and V in the limit fi — > with fiV held fixed is expressed 
by a matrix integral over a U(Nf) matrix Uq. The role of the QCD partition function is 
played by Z eS : 

Z cS (z) = f dU e^ TrUo , z = /iV± eS (8) 



The constant S e fj is a parameter inherited from the complete theory. This is often quoted 
in continuum contexts, but also holds on the lattice. It is in the latter sense that we write 
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the formula here. Under the conditions stated above, £ e g- = E. 

The quenched case is formally obtained by letting Nf approach zero. This additional limit 
is not very well defined since theories with a fractional Nj are non-local and therefore may 
develop additional singularities and needs special care 15| . The quenched case can also be 
studied by the "supersymmetric method", which appears more reliable |lo[[l^|. The upshot 
of such an analysis is that there are infrared renormalization effects from fluctuating fields 
that cannot be neglected. These effects can be absorbed into a redefinition of £ e ff ; which 
now becomes dependent on V } complicating the formulation of the crossover effect jl6| . 
Actually, E(V) diverges as V — > oo, but the divergence is probably weak enough to allow a 
separation of scales which maintains the applicability of the matrix integral formula for Z e &. 
It is somewhat unclear how solid this last claim is for finite N c . However, things simplify 
again if one also expands in The volume dependence of -^-£(V) enters then order by 
order in the jj- expansion. The leading term is volume independent, limN c ^ >00 ^-T l {V) = £ 
and the first subleading term has a logarithmic dependence on V. Thus, the dependence on 
V disappears also in the quenched case if we take the large N c limit. In short, the large N c 
limit - as we defined it - of ordinary quenched theories is free of infrared divergences and 
coincides with the large N c limit of theories with finite and fixed Nf > 1. 

The formula for Z e R would hold in any fundamental model that has the same flavor sym- 
metries and breaking pattern as QCD. However, the case of QCD enjoys one special property 
that we have not yet exploited, namely that fermions enter the action only bilinearly. This 
property has been a past source of other important results, like the geometrical structure 
of anomalies and mass inequalities. In our context, the bilinear structure of the fermionic 
action, when combined with chiral effective Lagrangian techniques leads to chiral random 
matrix theory. This produces a more QCD specific and a more detailed understanding of 
the crossover regime. 

The main result of chiral random matrix theory concerns the eigenvalues of the operator 
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171 ]. Because of the symmetry of the spectrum, and since we are working at zero 



topology, we can focus on A > 0, and order the eigenvalues as Ai < A2 < A3 < .... < \k, 
where K = 2N C L A . As the quantity N C L 2 increases and goes through a certain threshold 
more and more of the low eigenvalues attain universal distributions in terms of random 
variables z^, with Zk = A^SA^L 4 3^. The z^s are the square roots of the ordered eigenvalues 
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of a matrix C'C, where C is a complex random matrix of size R, distributed as a Gaussian: 

d^{C)=Me- RTrC ' c n d ^ (9) 

l<i<j<R 

At fixed lattice spacing the product N C L 2 can be increased by increasing iV c at fixed L. 
This produces the same effect as increasing L at fixed iV c . Once Ai and A2 are determined 
to be jointly distributed in accordance with random matrix theory, the scaling parameter 
EiV c L 4 can be numerically extracted. This shows in detail how large iV c reduction works 
in this case, since the parameter E is defined a priori in the limit in which N c and L are 
simultaneously infinite. 

Thus, our numerical task is as follows: We first need to show that it is enough to increase 
N c at fixed L to cause A12 to become distributed according to random matrix theory. Once 
this is established we should extract E for different values of the fixed parameter L, and 
show that the different numbers are the same, and therefore L- independent. Subsequently, 
we should check whether more traditional ways to extract E lead to consistent results. We 
have carried out the required numerical work and our results support the hypothesis of this 
paper. Details will be presented in section El 

Encouraged by these results, we now provide a somewhat different viewpoint. Until now 
we relied mainly on theoretical results that were obtained at finite N c . At finite N c it is 
unlikely that the entire Dirac operator is a random matrix in any precise sense. It is only 
a tiny fraction of its lowest eigenvalues that are distributed as if they were the eigenvalues 
of a random Dirac matrix. One could imagine framing the analysis as follows: Assume 
that one does the integral over all gauge fields subjected to constraints that fix TrA~ 2k for 

all K > k > 1. This seems to define a distribution of {Ai, A2, , Xk}- For finite N c the 

number K is a finite fraction of the total number of gluonic degrees of freedom. The totality 
of {Xk, K > k > 1} contains as much information as a full fledged (highly nonlocal) quantum 
field and therefore its fluctuations are not expected to have a structure largely independent 

of the detailed dynamics. It is therefore unlikely that the distribution of Ai, A2, , Xk have 

eigenvalue repulsion of the universal type for all adjacent eigenvalues. 

However, at infinite N c , K is a vanishing fraction of the total number of degrees of freedom 
since K is linear in N c and there are order N 2 gauge degrees of freedom. It is now more likely 
that the simplest possibility holds: as the size of A goes to infinity the repulsion between 
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all adjacent eigenvalues is of the universal type. At zero topology A has the structure 

/ C\ 

A =U o) (10) 

C is a square matrix of size K x K and depends parametrically on the fluctuating gauge 

(I \ 

fields. (We assume that 75 = ) . As such, the distribution of C inherits an exact 

VO -\) 

invariance under gauge transformations. The volume is finite, but this is not reflected by C 
being sparse in any precise sense. It is true that at infinite L entires of C that are far in 
site distance would be small: but this is not very relevant for finite L. Thus, the C-block of 
A is probably well described by a random matrix model with U (K) x U (K) invariance. By 
"well described" we mean that after "unfolding" (rescaling all eigenvalues by a local average 
level density) the spectrum of A will be universal in the limit N c — > 00. While the local 
average level density is gauge theory specific, bears all the dynamics, and in many respects 
is constrained by ordinary field theory universality, the local microscopic level distribution 
is unaffected by the gauge theory dynamics and obeys random matrix universality. We 
believe that many more matrix observables in large N c gauge theory exist, obeying such 
dual universalities. The Dirac operator is, we conjecture, a nice example of this concept. 
These conjectures lead to an effective weakening of the importance of space-time locality at 
infinite N c . 

It has been observed long ago that S^SB is equivalent to the statement that the density 
of eigenvalues of A is nonzero at very small eigenvalues, contrary to the free case, where it 
vanishes as |A| 3 (in four dimensions) At finite four- volume and finite N c this amounts 
to a roughly equally spaced spectrum near zero with eigenvalues and spacings of order jfj±- 
This kind of spectrum is obtained generically in RMT where eigenvalues repel and the 
spectrum tends to be rigid, with small fluctuations in the total number of eigenvalues in a 
given interval. Thus, RMT naturally provides S^SB. All one needs is that the randomness be 
sufficiently strong, and that the asymptotic order of magnitude of eigenvalues and spacings 
for large R, where R is the dimension of the the vector space the random matrix acts on, be 
given by Co ^ at - . This asymptotic behavior is very natural, since the basic unitary matrix V 
is unitary, with a spectrum restricted to the unit circle. This is what seems to happen in the 
I > l c phase of planar QCD on a torus of side I, where all large Wilson loops are expected 
to obey an area law. In short, at infinite N c , the occurrence of S%SB already at finite four 
volume is a plausible consequence of the randomness of the gauge fields. 
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The relevant formula from chiral RMT for our numerical work are the distributions of 
the two lowest eigenvalues and the distribution of their ratio in the Q = and Q — 1 
topological sector [l^]. The distribution of the lowest scaled eigenvalue Z\ and the second 
scaled eigenvalue z 2 in the Q = topological sector are given by 
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P2{z 2 ) = ]e~^z 2 r duu[ll{u) - h{u)h{u)\ (12) 
4 Jo 

The distribution of the ratio r = — in the Q = topological sector is given by 

P(r) = 4 (1 _ r2)2 I due W^)u 3 [I*(u) - h(u)I 3 (u)} (13) 

The distribution of the lowest scaled eigenvalue Z\ and the second scaled eigenvalue z 2 in 
the Q — 1 topological sector are given by 

1 _£ 

Pi(«i) = -e 4 ^i"^) (14) 

= V* / o 22 ^{^/ 2 (^)(^ 2 - *?)[i?(>/^-*?) - - 5) J ^( v^ 2 - *?)] 



4 Jo Zl 

+ 4h(z 2 )y/z$ - z\\h^zl - z\)h^zl - z\) - h^z\ - zl)I 2 ^zl - *?)] 

+ zll Q {z 2 )[ll^ zl - zl) - J X { % Izl - 2f)/ 8 (V 'zl - zl)]} (15) 
The distribution of the ratio r = — in the Q = 1 topological sector is given by 

I /'OC 2 , 

p(r) = — / d^ 3 e--{(l - r 2 )/ 2 (^)[/ 1 2 (2v / T 3 ^) - I (zVT^)I 2 (zVT^)] 
4r jo L 

+ Vl - H/i(z)[/ (zVl - r 2 )I 3 (zVl - r 2 ) - T/^zVl - r 2 )I 2 (zVl - r 2 )] 

+ / (^[/ 2 2 (^v / r^ :2 )-/i^v / r^ :2 )/3(2v / r^ :2 )]} (i6) 

The integrals involving modified Bessel functions were numerically evaluated to a high pre- 
cision and accurate plots of the various distributions were generated. 

IV. ESSENTIALS OF THE LATTICE FORMULATION. 

We used the simplest, single plaquette, pure gauge Wilson action, given by 

S= E Tr[U^{x) + Ul u {x)} (17) 
U^(x) = U,{x)U v {x + fj)Ul(x + u)Ut(x) (18) 
11 



We define b 



(3 

2M 



and take the large N c limit with b held fixed. As usual, b 



determines the lattice spacing a. The lattice is a symmetric torus of side L in lattice units. 
The gauge fields are periodic, x is a four component integer vector labeling the site, and fi 
either labels a direction or denotes a unit vector in the // direction. The link matrices U^{x) 
are in SU(N C ). The total topological charge, using the "overlap" definition |20j with mass 
parameter M set to —1.5 is denoted by the integer Q. 



Q = -Tr [H W {M)\ 



(19) 



We generated SU (N c ) gauge field configurations at several values of b. The evaluations of 
fermionic observables in a given gauge field background were separated by 25 updates. In 
our terminology, one update of the lattice corresponds to one Cabibbo-Marinari heatbath 
update for each one 

of the W- 1 ) SU(2) sub groups of SU(N C ) for each link, followed by 
one full SU(N C ) overrelaxation pass over the entire lattice [3|. Thermalization at a fixed L, 
N c and Q was achieved by performing a total of 500 updates. The starting configuration at 
Q — was 

f 1 for 1 < ii < L, 1 < i 2 , i 3 , «4 < L 



Z\ it = L, 1 < i 2 ,«3,^4 < L 

1 for 1 < i 2 < L, 1 < ii, i 3 , 24 < L 

z 2 i 2 = L, 1 < ii,z 3 ,Z4 < L 

1 for 1 < i 3 < L, 1 < ii, i 2 , 24 < L 

z 3 i 3 = L, 1 < ix,i 2 ,i4 < L 

1 for 1 < ^4 < L, 1 < ii, 22, ^3 < L 



(20) 



I z 4 z 4 = L, 1 < H,i 2 ,i 3 < L 
where z^, fi = 1,2,3,4 are randomly chosen members of Zn c . 

For a Q — 1 configuration we started with all links defined by the "uniform" instanton 
configuration below: 



up < 



e l : 

,27r(i 2 -l) 



e 

i; 

10; 

■ 27r(n-l) 

e : 



U, 



h3, 



H,l2,i3,U) 



.2^(^-1) 



e 

i; 
10; 



for i = j — 1; 1 < i 2 , i 3 , i i <L\ix = L 
for i = j = 3; 1 < 2 2 , 23, < L; ii = L 
for i = j ^ (1, 3); 1 < z' 2 , z' 3) «4 < £ 
elsewhere 

for i = j — 1; 1 < ii, i 2 , i 3 , L 
for i = j = 3; 1 < i 2 , i 3 , i±< L 
for i — j ^ (1, 3); 1 < i ls i 2 , i 3 , z 4 < L 
elsewhere 
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#3 J («i,«2,«3,«4) 



37r(i 4 -l) 
e L ; 



■ 2x(j 4 -l) 

e i 

i; 

10; 

. 27r(io-l) 

I — --4* — - 

e l 2 : 



21,22,23,24) 



e 

i; 
10; 



.2ir(io-l) 



(21) 



for i = j = 2; 1 < z'i, 22, 24 < L; 2*3 = L 
for 2 = j = 3; 1 < z'i, 22, 24 < L: 23 = L 
for 2 = j 7^ (2, 3); 1 < i x ,i 2 , 23, 24 < L 
elsewhere 

for z = j = 2; 1 < ii, z 2 , 23, 24 < L 
for z = j = 3; 1 < z'i, 22, 23, 24 < L 
for i = j ^ (2, 3); 1 < h,i 2 , 23, 2 4 < L 
elsewhere 

No tunneling during our thermalization cycle ever occurred, so the equilibrium configurations 
also have Q = 1. As is well known, during a simulation, the "overlap" topological charge will 
change only extremely rarely (far beyond what is observable in practice) whenever all gauge 
configurations produce a hermitian Wilson matrix H w that has a finite gap around zero. At 
large N c the gauge configurations come naturally out preserving this gap p and this explains 
the absence of tunneling events. In turn, the gap in the spectrum of H w is a consequence of 
a gap in the eigenvalue distribution of the parallel transporters round single plaquettes. The 
latter have a gap because we are working in a phase which is disconnected from the regime 
of lattice strong coupling by a phase transition where this gap forms. The existence of a gap 
in H w is of great practical value when dealing with the overlap Dirac operator because it 
eliminates the need to deal separately with the subspace of H w associated with very small 
eigenvalues in absolute magnitude. Technically, we do not need to "project out" these states 
when we use the pole approximation to evaluate the action of the unitary matrix V j^j]. 
Also, a robust gap in H w makes the definition of Q practically unambiguous (in the sense 
that one can vary the parameter M in a reasonable range and the assignment of topological 
charge will not change). Note that at large N c the computer cost balance between the single 
pass j^lj and double pass [2^ switches relative to N c = 3 and, typically, the single 
pass version, with Zolotarev coefficients J^, is more efficient. Single pass with too large N c 
would be prohibitive because of memory considerations, but, in practice, we did not need 
to deal with this issue. 
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TABLE I: List showing the number of configurations used at different L, N c and Q for the analysis. 
The table also displays results for the average of the ratio of the first eigenvalues to the second 
eigenvalue divided by its RMT value and the estimates of the condensate from the first and second 
eigenvalues. 



b 


L 


N c 


Q 


No. of conf. 


(Ai/A 2 ) 
(A1/A2/RMT 


1 




0.346 


9 


11 


l 


96 


0.99(3) 


0.1551(20) 


0.1553(10) 


0.350 


6 


13 





436 


1.99(2) 


0.1065(05) 


0.1356(04) 


0.350 


6 


17 





511 


1.65(2) 


0.1108(08) 


0.1333(05) 


0.350 


6 


23 





446 


1.11(2) 


0.1421(12) 


0.1485(06) 


0.350 


6 


29 





599 


1.06(2) 


0.1419(10) 


0.1450(05) 


0.350 


6 


37 





287 


1.01(3) 


0.1458(15) 


0.1464(07) 


0.350 


6 


37 


1 


192 


1.05(2) 


0.1440(13) 


0.1469(07) 


0.350 


6 


43 





292 


1.02(3) 


0.1401(14) 


0.1406(07) 


0.350 


7 


17 





346 


1.07(3) 


0.1376(13) 


0.1412(06) 


0.350 


7 


19 





315 


0.99(3) 


0.1407(14) 


0.1408(06) 


0.350 


7 


23 





440 


1.01(3) 


0.1415(12) 


0.1426(05) 


0.350 


7 


29 





348 


1.02(3) 


0.1441(13) 


0.1454(06) 


0.350 


7 


29 


1 


288 


1.04(2) 


0.1403(09) 


0.1426(05) 


0.350 


8 


13 





310 


1.04(3) 


0.1352(13) 


0.1371(06) 


0.350 


8 


17 





270 


1.04(3) 


0.1384(15) 


0.1399(07) 


0.350 


8 


23 





257 


1.03(3) 


0.1398(14) 


0.1413(07) 


0.350 


8 


23 


1 


288 


1.01(2) 


0.1418(10) 


0.1426(06) 


0.350 


10 


11 





64 


1.05(6) 


0.1313(27) 


0.1333(17) 


0.355 


8 


23 





288 


1.03(3) 


0.1192(13) 


0.1205(05) 


0.355 


9 


17 





288 


1.04(3) 


0.1204(12) 


0.1220(06) 


0.355 


10 


13 





288 


1.09(3) 


0.1118(11) 


0.1160(05) 


0.3585 


9 


17 





336 


1.10(3) 


0.1083(10) 


0.1121(04) 


0.3585 


9 


23 





300 


1.03(3) 


0.1062(11) 


0.1078(05) 
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V. OVERVIEW OF NUMERICAL WORK 



We started our study at a lattice gauge coupling b = 0.350 where the bare 't Hooft 
coupling is given by fS. = \. At this lattice coupling a lattice size of L - 6 is very 
close to critical Thus, a, the lattice spacing, is approximately given by a = l f. To get 
some physical idea of this size in QCD terms consider N c = 3 and assume l c = f-, where 
T c is the pure gauge finite temperature gauge transition. This means that in QCD terms 
a = g^r- ~ Q-75 fermi _ | j erm ^ We kept ^ fixed at this value and simulated the pure gauge 
sector at L — 6, 7, 8. From the continuum viewpoint this provides for increasing physical 
volumes at fixed physical lattice spacing. 

At each L we simulated several increasing values of N c , until we reached a regime where 
the calculated Ai < A2 had a joint RMT distribution. This question was addressed in a scale 
independent way by considering the distribution of the ratio r = p{r)dr. The random 
variable r takes values in the segment (0, 1) and the behavior of the distribution near 1 
reflects the typical RMT eigenvalue repulsion. The behavior of the distribution near is 
dominated by the universal features of "edge" behavior in RMT. 

The point of looking at p(r) is that the defining property of randomness in the context 
of random matrix theory is eigenvalue repulsion. The tail of the distribution of p(r) near 
r = 1 reflects this repulsion best. As N c increases at fixed L, the lowest eigenvalues of 
—A 2 enter the regime where they ought to be described by chiral random matrix theory if 
there is S^SB. They enter the regime by being gradually squeezed closer together by the 
increase in N c (recall that the matrix V is unitary, so its eigenvalues are restricted to the 
unit circle), until their distribution is totally governed by eigenvalue repulsion. This means 
that we expect that the approach to RMT will be by distributions p(r) which show a gradual 
depletion close to r = 1, approaching the RMT prediction from above in that region, and 
therefore, to compensate for the fixed normalization, from below in the region of small r 
close to zero. This means that the average of r will approach its RMT limit from above as 
N c increases. 

Therefore one expects that once the empirical average of r is close to its RMT limit, the 
entire distribution p(r) will fall into statistical agreement with the RMT prediction. This 
expectation was confirmed by our data. Note that using the random variable - instead of 
r is ill advised, since then the behavior at r = gets emphasized, and the variance of - 
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diverges. We found that we needed at least several hundred gauge field configurations to 
be able to make a reasonable determination that the RMT regime has been entered. In 
general, the fluctuations are quite large and our statistics are relatively low. To be safe, we 
feel one should multiply our errors by a factor of 2 to 3 in order to get some reasonably 
reliable confidence intervals. To check whether we are not missing some hidden effect here, 
we went to non-zero topology and checked whether the method of random matrix theory 
works similarly and consistently there. The results at nonzero topology indeed behaved as 
expected. 

To this point we have not introduced a mass yet. To see how the large N c limit sets in 
numerically we looked for a quantity that approaches S as fi — > 0. The simplest such quantity 
is the bilinear condensate but, for nonzero \i this quantity diverges quadratically 

when a goes to zero. This effect indicates that we may need to subtract some large numbers 
from the measured value, enhancing the noise/signal ratio. We preferred a quantity with 
a milder behavior in the ultraviolet. We chose the zero momentum pion-pion scattering 
matrix element. At nonzero but small quark mass and at leading order in the expansion 
it is given by the gauge average of Tr ( — rj- 2 ) . One has then: 




While the quantity we are looking at would still have a divergent contribution as a goes 
to zero, that divergence is only logarithmic, and not threatening numerically. Moreover, 
the numerical effort involved in evaluating the latter quantity is equal to that involved in 
computing the former. We found consistent results, but also learned that the random matrix 
method is by far more accurate as a numerical tool to establish both the presence of S^SB 
and the value of the condensate. 

The comparison between the direct determination of the condensate and the one via 
random matrix theory left us somewhat uneasy because the direct determination, in itself, 
did not provide convincing evidence for S^SB. To be sure, if one first postulated S^SB, 
the numbers were consistent with the findings using random matrix theory. Nevertheless, 
to reassure ourselves that nothing was wrong, we looked at identical issues in the 't Hooft 
model. (Although the model is in two dimensions, infinite N c brings it closer to the planar 
four dimensional case, because spontaneous symmetry breaking of continuous symmetries is 
now an open option.) In the 't Hooft model we can calculate analytically the value of the 





(22) 
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condensate in the continuum limit. This provides an independent check showing that our 
methods do yield correct results and indeed one can trust the random matrix method even 
when direct measurements might appear inconclusive by themselves. 

We continued the zero mass analysis to finer lattices, looking at couplings b = 
0.355, 0.3585. These couplings are chosen to be almost equal to the critical bare coupling 
b c (L) for sizes L = 7, 8 3]. Thus, one expects the lattice quantity Es(6 c (L)) to go as j so 
long as the variation of the normalization constant for the pseudoscalar density varies slowly 
over the range. 

We can estimate from 1-loop continuum perturbation theory how much one expects the 
renormalization constant to contribute to the variation of with lattice spacing. To get 
a finite continuum condensate in some renormalization scheme one needs to multiply the 
lattice condensate by an appropriate constant Zg{X)- Slightly modifying known results, we 
have to 1-loop: 

Z S (X) = 1 + \[-^ log(a//) + c] (23) 

Here, /x is a renormalization point and a is the lattice spacing, c is a scheme dependent 
constant 1271, irrelevant for the estimates below. Thus, 



dZ, 



da 



(24) 



For us the lattice spacing goes as -j-, therefore ^ ~ 1, and b ~ 0.355 results in a 5Z S ~ 0.008. 
The contribution coming from 5 A also is small. The variation SZ S is away from a value of 
order unity, so contributes to the change in the condensate at the level of one percent. 
But, the condensate's non-anomalous scale dependence is a 3 which would amount to a non- 
anomalous change of about 60 percent. Within our accuracy the anomalous contribution 
coming from Z$ is therefore too small to be detectable. Thus, scaling will be verified just 
by observing the lattice go approximately as 

VI. NUMERICAL RESULTS IN FOUR DIMENSIONS 
A. Analysis based on chiral random matrix theory 

The two lowest non-zero eigenvalues of interest of the massless Dirac operator were com- 
puted using the Ritz functional method 2^ on D Q D\. We worked in a fixed chiral sector 
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and computed the two lowest eigenvalues of D D\. When Q = 1 we made sure to work in 
the chiral sector where there is no zero mode. 

The two lowest eigenvalues were used to obtain an estimate of the chiral condensate in 
the large N c limit. Let Ai < A2 be the two lowest eigenvalues in a fixed background gauge 
configuration. For large enough N c the probability distribution of the ratio r = ^ is given by 
a universal function that only depends on Q and is defined by chiral random matrix theory. 
If the probability distribution of the ratio agrees with the prediction from the chiral random 
matrix theory, one is confident that the distributions of the scaled eigenvalues Z{ = AjEjiV c L 4 
should also agree with the predictions from chiral random matrix theory. If this confidence 
is well based we must find that E x equals E 2 . 

For a given L, N c and Q we defined E; as the numbers that made the observed average 
of Zi agree with the numbers from chiral random matrix theory with a Gaussian weight as 
explained before. We checked to see if Ei = E 2 and if the detailed scaled distribution agreed 
with chiral random matrix theory predictions. The list of different values of L, N c , b and 
Q where we did simulations is shown in table HI For most cases, we see general agreement 
with chiral random matrix theory. 

Fig. ^ compares our data to chiral random matrix theory for the distribution p(r) in the 
Q = sector for lattice size 6 4 and coupling b = 0.350 as N c increases. One should keep in 
mind that what we are plotting are histograms, a collection of integers obtained by binning 
the data. The integer associated with a given bin undergoes Poissonian fluctuations, so the 
standard deviation is equal to the square root of the same integer. For this reason, the more 
populated a bin is the larger is the absolute error bar associated with it. The distributions 
in the different bins are somewhat correlated, but there are sizable fluctuations from bin to 
bin and this is expected. Even within our limited statistics it is obvious that a dramatic 
change in the distribution of the ratio occurs between N c = 17 and N c = 23, and that by 
N c = 29 the distribution has stabilized, essentially reproducing the infinite N c limit. 

We now look at the distributions of the rescaled two eigenvalues as a function of N c , still 
at L = 6 and b = 0.350 in Fig. Eland in Fig. El We see that merely looking at one eigenvalue 
in isolation may be misleading (compare N c = 17 to N c = 23 for instance). A glance at 
table U shows that the determinations of Ej differ significantly between i = 1 and i = 2 
for N c = 13, 17 at L — 6, b = 0.350, confirming the indication form the ratio distribution. 
Again, the scatter in the data is consistent with Poisson statistics. 
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We now pursue the question of L dependence by keeping b = 0.350 and looking for the 
infinite iV c limit of the distributions and correctly rescaled eigenvalues for higher values 
of L. In addition, we wish to ascertain that these distributions depend on topology in 
the precise way RMT dictates. Moreover, as table indicates, the numerical values of the 
Sj are topology independent within our statistics. The results we would like to look at 
at this point are collected in Fig. |U (for the ratio distributions), Fig. El (for the smallest 
eigenvalue) and Fig. El (for the next smallest eigenvalue). Our Q — 1 simulations focused 
on a value of N c large enough that the Q = case is consistent with RMT already. In a 
marginal case consistency with RMT may not hold for Q = 1 even if it holds for Q — 0, 
because the two nonzero lowest eigenvalues at Q = 1 are repelled by the zero eigenvalue and 
therefore somewhat higher in magnitude than in the Q = case. The statistical features 
and fluctuations at Q — 1 are similar to the ones at Q — 0. 

We now pursue other values of b, at appropriate volumes and N c values where RMT holds 
and we can focus on the numerical value of E, ultimately checking consistency with asymp- 
totic freedom scaling. A sample of results is collected in Fig. (for the ratio distributions), 
Fig. El (for the smallest eigenvalue) and Fig. El (for the next smallest eigenvalue). 

We collect our determinations of the chiral condensate at b = 0.350 in Fig. El The data 
points displayed in this figure show that, within errors, Si = E2. The errors we displayed 
for the Ej are obtained by assuming that the law of large numbers holds and the average 
of the individual eigenvalue is Gaussianly distributed. Our statistics are somewhat meager 
and therefore it is safer to view our numbers as reliable to within two or three standard 
deviations. 

To get some feeling for the physical order of magnitudes we are getting we consider the 
y-axis in Fig. we see that E^ has a range of 0.142 ± 0.006. Ignoring the scalar density 
renormalization factor we can get an approximate number for the continuum value of the 
condensate. 

rs(6 = o.35)l 13 _ oog8±0 004 (25) 



2\M\ 

where M is the Wilson mass that enters the overlap Dirac operator. We have set M = —1.5 

in our calculations. Taking into account that the lattice critical size L c (b) ~ 6 at b = 0.350 

and assuming that l c T c = 1, we get E^ ont (6 = 0.35) = 0.588T C . Taking T c « O.QtJcf « 
I 1/3 

26J gives S^ cont ~ 155MeV. If we now assume that N c = 3 is large enough to 



264MeV 
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ignore ^- corrections, {ipip) n c =3 ~ (224MeV) 3 . In spite of the roughness of our assumptions, 
the number comes out close to experiment, but one should not make too much of this, 
because the renormalization factor can easily change the scale by as much as 25%. 

One can try to get a rough feel for the latter effect by using tadpole improved one loop 
estimates J3] for Zs in the MS scheme 

ZgSQi) = 1.40{1 - M^[61ogM - 1.33]} (26) 

The prefactor 1.40 is \^ d , where the tadpole-improved [2^ Wilson mass is given by 

_ M ttd = |M| tad = 1^1 -^-"0) (27) 

u 

uo is given by the fourth root of the average plaquette, a number quite close to unity. The 
entire one loop correction is numerically negligible for /i ~ 2 GeV, so Z| 1S (2 GeV) = 1.40. 
This increases the scale by about 12%, to ^-(ipip) « (174MeV) 3 . For QCD this trans- 
lates into (V^qcd ps (251MeV) 3 , close to results obtained in quenched lattice QCD [30]. 
Obviously, before a scale can be reliably quoted one would need to determine Z s outside 
perturbation theory, by more numerical work. 



B. Extracting the chiral condensate from mass dependence 

The conventional method to estimate the chiral condensate is to obtain a stochastic 
estimate of the fermion bilinear (i/ji/j) (/x) as a function of fermion mass \x and perform an 
extrapolation to zero /i. Since, strictly speaking, we define the large N c limit of the massless 
theory by first keeping fi finite while taking N c to infinity, and only subsequently taking 
fi to zero, the conventional method is the most conservative way to go about establishing 
S^SB and estimate the order parameter ( , 0'0}(O). However, it is not practical to make N c 
large enough to disentangle the effect we are after from finite N c effects. As we shall see, the 
RMT route is much better in practice at determining the infinite N c behavior from finite iV c 
results. 

We evaluate (b\A^ 1 ( / u)yl _1 (/x)|6) for a chiral "source" vector \b >, randomly drawn from 
a Gaussian distribution, and obtain a stochastic estimate of the trace we are interested in: 

F 1 ( l x) = TrA^\ f i)A- 1 ( f x) (28) 
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We then estimate the condensate from 



E^KmFrCu). (29) 

Using one random source per configuration, we computed Fi(fi) on 48 different thermal- 
ized gauge background configurations at three different values of (L, N c ). These simulations 
were all done at our chosen coupling, b = 0.350. The results are plotted in the three left 
panels of Fig. ^2 and compared with the estimate of E 1 ' 3 from the previous subsection. 
The stochastic estimate seems to be consistent with the estimate from chiral random matrix 
theory but the range covered on the y-axis is too large to make this claim convincing. 

A more careful examination reveals that much of Fi(fi) comes from a free field contribu- 
tion and one could subtract that free field contribution by an explicit computation, without 
affecting the /i — > limit so long as there is S^SB. The coupling is not weak, so one must 
use the "right" free field case: The natural choice is to employ the tadpole-improved Wilson 
mass in H w , when computing the free fermion contribution. Denoting the free field result 
by F{(fi) we carried out the subtraction and the result is plotted in the three right panels 
of Fig. El We see that indeed Fx (fi) is dominated by the free fermion contribution and that 
the subtracted quantity still is consistent with the chiral random matrix theory prediction 
for S 1 / 3 . 

As explained in Section |Vj one can improve the situation by defining a quantity that has 
less of a divergence in the ultraviolet: 

F 2 {ii) = 2fx 3 Tr[A^\fx)A-\fi)f (30) 

It is easy to see that 

limF 2 1/3 (/,) = E 1 /3 (31) 

This is a particular attractive alternative, since a stochastic estimate of -^(/-O needs no more 
numerical computation than what is needed for Fi(/j,). We note that 

*M"*£P?i) (32) 

and therefore the ultraviolet divergent term, which is linear in /i and additively contributes 
to Fi(fi), does not contribute to .F2 (/•*)• The first sub- leading term in F 2 (fi) is proportional 
to /i 3 and therefore we expect a better estimate of S 1//3 from F 2 (fx). The results for F 2 (fx) 
are plotted in the three left panels of Fig. [T21 these results are obtained simultaneously 
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with those for at no additional numerical cost. Again, the results seem consistent 

with the prediction obtained using chiral random matrix theory in the previous subsection. 
Similarly to Fi(fi), there is a "free field" contribution to F 2 ([i) and we subtracted this "free 
part", F 2 (fi) (with tadpole improved parameters), by a direct computation. The resulting 
subtracted quantity is plotted in the three right panels of Fig. H21 and one can see that this 
time the subtraction had no discernible effect. As expected, the subtraction of the free field 
expression mainly affects the ultraviolet contribution. Consistency with the RMT results is 
better for F 2 (fi) than for but the RMT approach still appears significantly superior 

in practice. 



C. Approximate scaling 

As explained, the coarsest evidence for scaling should come by simply observing that the 
quantity s(b) = L c (b) (Jnm) 3 stays constant as b is varied. The lattice critical size L c (b) has 
been determined numerically in [^] to be given by: 

, 11 x 51/121 

L c (b) = (0.260 + 0.015) —^—) e 24 ^ 11 



48tt 2 6 7 , 
,& 2 - 0.589646 + 0.08467 

1 = b 2 - 0.502276+ 0.05479 ^ ' 

This formula uses tadpole improvement and provides an interpolation valid for 6 < L c (b) < 
10. For a given b the relevant calculations are carried out at L > L c (b). 

Looking at the table we find s(0.350) = 0.59(6), s(0.355) = 0.58(6) and s(0.3585) = 
0.58(6). Even a limited investigation on a very coarse lattice gives a consistent number: 
s(0.346) = 0.55(5). If we include the factor Z^{2 GeV) = 1.40, we obtain a continuum 
statement at infinite N c : 

^r(W) m (2 GeV) w (0.65) 3 (34) 
VII. NUMERICAL RESULTS IN TWO DIMENSIONS 



Two dimensional large N c QCD |31[ can be used as a soluble example for the phenomenon 
of S^SB at infinite N c . In the absence of exact zero modes, just as in four dimensions, we 
expect RMT to work; moreover, it is the same RMT model that applies. Also, one can define 
in entirely analogous ways the functions Fi,F 2 . The major advantages are: simulations are 
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relatively very fast, enabling one to compute the entire spectrum of A and, in addition, the 
continuum value of £ is known from the analytical solution of 't Hooft. 

For completeness, we start with a brief derivation of the exact value of the chiral conden- 
sate in large N c QCD P, [y| . Taking care of the wavefunction normalization, we have 

2 

TrA^\fx)A-\fi) =4M 2 ^| (35) 

Mo 

in the limit of small quark mass, where only the first term in the infinite sum over poles is 
expected to contribute (for this to be true the infinite sum over poles needs to be regulated 
to eliminate ultraviolet divergences - because the underlying model is renormalizable, it does 
not matter how exactly this regularization is done) and ro is the corresponding residue, /ig 
is the "pion" mass squared. The quark mass is related to \i by m q = 2|M|/i. The quark 
mass in units of the gauge coupling is given by a dimensionless parameter 7: 



Using (J35|) and (J36|) . we get 



V2^b 

2|M| v / 7 > o 2 
\f2ixb \i\ 

The lowest pseudoscalar mass is given by 

.2 _ 2 ^Vi 



™ q = -5(2= (36) 



(37) 



Mo = (38) 



and its residue is given by 

2 71 
r °=3 

Putting together all the information, we get 

^ \M\ 



(39) 



(40) 



We will use M = — 1 and b = 1 for our numerical analysis. Previous experience Jfj] tells us 
that b = 1 is sufficiently large to expect the effects of finite lattice spacing to be relatively 
small, an expectation confirmed by our findings here too. This provides a numerical value 
for the lattice condensate: X = -A= = 0.2303. 



The numerical computation using overlap fermions in two dimensions was done in a 
manner similar to the one in four dimensions. However, we performed a full diagonaliza- 
tion of the overlap Dirac operator in two dimensions, eliminating this potential source of 
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statistical noise (in practice, the fluctuations associated with the stochastic nature of the 
trace evaluation are substantially smaller than those associated with varying the gauge field 
background). Thus, we did not resort to a stochastic estimate of Fi(fi) and F 2 (/j,) in two 
dimensions and, as we shall see presently, obtained data very similar to four dimensions, 
where we are forced to use a stochastic estimator for the trace. We generated a total of 1024 
configurations at a fixed L and N c . We worked with the following (L,N C ) pairs: (4,57), 
(5,47), (6,37) and (7,31). All configurations were gauge field backgrounds with no exact 
fermion zero modes. 

Like in four dimensions, we start with a plot of the distribution of the ratio of the first 
eigenvalue to the second eigenvalue and compare it with chiral RMT. This is shown in the 
four left panels in Fig. ED As in four dimensions, the agreement is fairly good and one sees 
large scatter at values of r where p(r) is large. The four middle panels and the four right 
panels in Fig. El show the comparison between data and chiral RMT for the two lowest 
eigenvalues. In relating the chiral RMT variable z to the eigenvalue A of the overlap Dirac 
operator via z = AEiV c L 2 , we used the known result for the condensate, namely, E = 0.2303. 
Our main point is that we know in this case that S^SB occurs and the plots still have the 
same general structure as in four dimensions where we wanted to determine whether S^SB 
takes place. 

Since we have computed all the eigenvalues of the overlap Dirac operator, we can also 
obtain a plot of the spectral density, p(A). This is impossible to do in four dimensions, except, 
maybe, for small values of N c . The result is in Fig. EJ where we show the distribution only 
close to A = 0. Getting the eigenvalue density even in this limited range would be very costly 
in four dimensions. As is well known p(0) = H/tt and we see that the theoretically 
expected number, 0.0733, is correctly reproduced, with an accuracy and consistency similar 
to that obtained in four dimensions, using RMT. 

We can try to see if one can also get the chiral condensate from a direct extrapolation of 
Fi(fjL) and F 2 ([i). Plots of Fx(p) and F 2 (p) along with their corresponding result after the 
subtraction of the tadpole improved free field contribution are shown in the various panels 
of Fig. EJand Fig. El We see that these more direct observables indeed are poorer practical 
indicators for S^SB and give less accurate estimates of the condensate if S^SB is assumed. 
The overall structure of the data is similar to four dimensions, indicating that nothing is 
wrong with our interpretation of the four dimensional results. 
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VIII. SUMMARY 



Our evidence supports the basic hypothesis of this paper: that, in the planar limit, 
continuum SU(N C ) gauge theory, defined on an Euclidean four dimensional torus of side I, 
breaks chiral symmetry spontaneously so long as I > l c , where l c is a physical length of order 
1 fermi in QCD units. 

Some may find it surprising that a continuous symmetry is spontaneously broken in a 
finite volume. Maybe the following rather trivial observations would convince the skeptics 
that this is not that unexpected. Suppose we consider massless QCD, in a Hamiltonian 
formulation, at zero temperature, defined on S 3 , where the scale of the S* 3 is s. It is easy 
to see that, as s — > 0, ('ip(x)ip(x)) ~ with a nonzero constant. The much harder 
question is what happens as s — > oo: perturbatively the condensate vanishes, but, if there is 
S^SB, it should drop to a nonzero constant and level off. Now, make the temperature finite. 
This causes the condensate to disappear, no matter how small the non-zero temperature is. 
However, if we first go to the planar limit, the condensate will not disappear at finite and 
small temperatures. 

Our investigation has led us to a picture of spontaneous breaking of chiral symmetry 
at infinite N c which makes the phenomenon appear generic. In other words, one almost 
ends up with the conclusion that the more difficult problems are to explain how, in more 
complicated cases, chiral symmetry fails to break spontaneously in the planar limit. 
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FIG. 1: The distribution of the ratio of the first eigenvalue to the second eigenvalue in the zero 
topological sector for various N c at L = 6 and b = 0.350 is compared with the prediction from 
chiral RMT in four dimensions. 
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FIG. 2: The distribution of the first eigenvalue scaled by the corresponding average value in the 
zero topological sector for various iV c at L = 6 and b = 0.350 is compared with the prediction from 
chiral RMT in four dimensions. 
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FIG. 3: The distribution of the second eigenvalue scaled by the corresponding average value in 
the zero topological sector for various N c at L = 6 and b = 0.350 is compared with the prediction 
from chiral RMT in four dimensions. 
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FIG. 11: The three panels on the left side show the stochastic estimate of F l (/j.) in four di- 
mensions for three different combinations of L and N c . The three panels on the right side show 
the same estimate after the subtraction of the tadpole corrected free field contribution. The solid 
circle with errorbars is the estimate for the chiral condensate from chiral RMT. 
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FIG. 12: The three panels on the left side show the stochastic estimate of F 2 (/u) in four di- 
mensions for three different combinations of L and N c . The three panels on the right side show 
the same estimate after the subtraction of the tadpole corrected free field contribution. The solid 
circle with errorbars is the estimate for the chiral condensate from chiral RMT. 
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FIG. 16: The four panels on the left along with the four panels on the right show plots of i^O^) 
in two dimensions before and after the subtraction of the free field contribution. The solid square 
present in all the panels is the exact continuum result for the chiral condensate. 
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